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Within the framework of the instantaneous Bethe-Salpeter equation, we present a detailed analysis 
of light meson spectra with respect to various parameterizations of confinement in Dirac space. 
Assuming a linearly rising quark-antiquark potential, we investigate two different spinorial forms 
(Dirac structures), namely |(1®1— 7°<8>7°) as well as the ?7./i(l)-invariant combination §(I® 1 — 7 5 ® 
7 5 — 7 M ® 7 M ), both providing a good description of the ground state Regge trajectories up to highest 
observed angular momenta. Whereas the first structure is slightly prefered concerning numerous 
meson decay properties (see Q]), we find the [/A(l)-invariant force to be much more appropriate 
for the description of a multitude of higher mass resonances discovered in the data of the Crystal 
Barrel collaboration during the last few years. Furthermore, this confinement structure has the 
remarkable feature to yield a linear dependence of masses on their radial excitation number. For 
many experimental resonances such a trajectory-like behaviour was observed by Anisovich et al. 
We can confirm that almost the same slope occurs for all trajectories. Adding the [/A(l)-breaking 
instanton induced 't Hooft interaction we can compute the pseudoscalar mass splittings with both 
Dirac structures and for the scalar mesons a natural mechanism of flavour mixing is achieved. In 
the scalar sector, the two models provide completely different ground state and excitation masses, 
thus leading to different assignments of possible qq states in this region. The scalar meson masses 
calculated with the structure |(I ® 1 — 7 ® 7 s — 7 M ® 7 M ) are in excellent agreement with the 
K-matrix poles deduced from experiment by Anisovich and coworkers. 



I. INTRODUCTION 



In the last few years, evaluations of the Crystal Barrel pA^-annihilation data set provided a lot of new meson 
resonances in the mass region 1000 - 2400 MeV fj]||[§§|. The experiment was performed with the Crystal Barrel 
detector at LEAR where antiproton-proton annihilation into Tr°nr), ir°ir°ri, 4tt°, Tr°n and Tr°n' has been studied 
up to an incident beam momentum of 1.94 GeV/c. 

From a theoretical point of view, higher excitation resonances in the meson mass spectrum are of great interest 
because they reflect the underlying confinement structure at the quark level. However, from first QCD principles 
especially the spinorial Dirac structure of the confinement is largely unknown. Therefore a phenomenological study 
of different Dirac structures in a fully relativistic framework seems worthwhile. From earlier works ]^,|l0| it is known 
that a linearly rising confinement potential with a Dirac structure of the form 1 ® 1 — 7 ® 7 provides a very good 
description of the Regge trajectories up to total angular momentum J = 6. However, the suitability of a confinement 
parameterization should not only be linked to the bulk of experimental ground state masses in each channel but also 
to their corresponding radial excitations. For that, a multitude of well established higher resonances has to be known. 
Unfortunately this is not the case so far, although especially during the last two years the situation has improved. For 
instance, while the Particle Data Group 98 (PDG 98 JO]) did not list any radial excitations for the isovector states 
ai(ll ++ ), 02(12++), 04(14++) and &i(lU ), in their latest issue (PDG 00 |47|]) they state three new such resonances: 
the vector meson ai(1640) observed by the authors i n M ( a similar resonance, the <2i(1700), is stated in ||), and two 
tensor mesons, the 02 (1660) H|| and the 02 (1750) | ]48|j4g| ] . Furthermore a lot of resonances have been found in the 
last few years, which do not appear in the latest PDG-listing. For example: o (2025) @, oi(2100) |J, oi(2340) [|, 
o 2 (2100) @, a 4 (2260) @, ^(2040) §, ^(2300) §, A(1700) @, /i(2340) §, / 4 (2320) §. We compare these new 
experimental resonance positions with the eigenvalue spectrum calculated with the Dirac structure 1 (g> 1 — 7 <£> 7 
(model A) on the one hand and the structure 1 <8> 1 — 7 5 <8> 7 5 — 7 M <8> 7^ (model £>) on the other hand. Whereas 
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the first structure produces excitation masses which overestimate the corresponding experimental masses up to 400 
MeV, the second structure yields a remarkably good agreement with the newly observed resonances. Here, deviations 
are in general less than 100 MeV. Furthermore, in model B the squared masses of the resonances show, in contrast 
to model A, a linear dependence on their radial excitation number, M 2 cx n, very similar to the behaviour of many 
experimental resonances recently observed by A. V. Anisovich, V. V. Anisovich and A. V. Sarantsev M. 

It turned out that the comparatively strong coupling of positive and negative energy components of the Salpeter 
amplitudes in model B is responsible for this desired lowering of the excited meson masses. We will illustrate this 
effect by carrying out the nonrelativistic reduction of the full Salpeter equation for both models. Whereas even in 
this limit the above mentioned coupling does not vanish for model B, this reduction provides the usual Schrodinger 
equation for model A, i.e. the complete decoupling of positive and negative energy components. Indeed, model A 
reduces in this limit to a particular version of the nonrelativistic quark model (NRQM) that provides a satisfactory 
description of both meson and baryon ground state masses 0]. Of course, due to the absence of the negative energy 
components, the NRQM in [Q fails for higher resonances as does the often cited model of Godfrey and Isgur 

The relativistic Salpeter framework presented here contains a priori all small desired spin-orbit terms in order to 
describe the small mass splittings between states that can be attributed to the same orbital angular momentum 
multiplets. Examples are oi(1260) and a 2 (1320), /i(1285) and / 2 (1270) and A' 2 (1770) and A 3 *(1780). One should 
emphasize that the low positions of the a (980) and /o(980) cannot be explained by these intrinsic spin-orbit terms 
alone. In order to describe these splittings an additional residual interaction has to be adopted. We will use an 
instanton induced effective quark interaction discovered by 't Hooft |L6|], which also accounts for the correct vector- 
pseudoscalar and K-ir-rj-T]' mass splitting. Whereas in the nonrelativistic model ||l4j| this interaction only acts for 
pseudoscalar mesons, in the fully relativistic Salpeter model jl7],[L8| it also acts in the scalar sector and provides a 
possible interpretation of the still unknown scalar ground state nonet [ [l9|| . However, as we will see this interpretation 
strongly depends on details of the confinement force, in particular on whether the used Dirac structure induces 
additional spin-orbit forces or not. In fact, as the Dirac structure 1 (g> 1 — j 5 <E) 7 5 — 7 M <E> "f/j, generates such additional 
spin-orbit terms, model B provides a completely different result for the scalar mass spectrum compared to the earlier 
computations in model A which produced an almost flavour singlet /o-state at approximately 1 GeV and the flavour 
octet states /o, Aq, ao in a mass region around 1.4 GeV (see |l9|): In model B, the calculated scalar ground state 
masses (and also their radial excitations) are roughly 200-300 MeV lighter. Moreover . they show a remarkable 



agreement with a lot of K-matrix-poles deduced by V. V. Anisovich and coworkers [|22| g4 43 45|] from experiment. 

This paper is organized as follows: In section II we briefly comment on the Bethe-Salpeter equation for a quark- 
antiquark bound state and display the approximations (instantaneous approximation, free quark propagators) that 
lead to the Salpeter equation, which constitutes the basic equation of our model. Section III is devoted to an extensive 
discussion of the effects of the model interactions adopted on the description of the experimental mass spectrum: In 
part A, focusing on the confinement force alone, we compare the resulting bound state masses of model A and model 
B to the complete J > mass spectrum and especially to the recently observed higher resonances quoted above. 
The nonrelativistic reduction of both models is presented in part B of section III, where we also define the positive 
and negative energy components of the Salpeter amplitude. A brief discussion concerning spin-orbit effects in both 
models is given in part C. In the last part of section III, we focus on the pseudoscalar and scalar mesons by adding 
the residual 't Hooft interaction which only acts for mesons with vanishing total angular momentum. There we will 
also comment on a set of K-matrix poles found by V. V. Anisovich and coworkers [p2]- p4] , fl^ |4"5|] by comparing to our 
results. Finally we give a summary and conclusion in section IV. 



II. A COVARIANT QUARK MODEL IN THE INSTANTANEOUS BETHE-SALPETER APPROACH 



In quantum field theory a quark-antiquark bound state with four-momentum P and mass M , M 2 = P 2 , is described 
by the Bethe-Salpeter (BS) equation for two fermions fl2| . In momentum space, this equation reads: 

X P (P) = S[ (| +p) J 0-; [- lK {P^p') X p {p')\ SlVf +P), (1) 

where p is the relative four-momentum between the quark and the antiquark, K denotes the infinite sum of their 
irreducible interactions and the corresponding full Feynman propagators are labeled by S[ and S^, respectively. The 
BS amplitude \ P is defined in coordinate space as the time-ordered product of the quark and the antiquark field 
operator between the bound state \P) and the vacuum: 
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Xap(xi,x 2 ) :=(0 \T ^{xx^Hx^ P) = e 



-iP-(a;i+X2)/2 



d 4 p 

(2^ 



(2) 



where a and /3 are multi-indices for the Dirac, flavour and colour degrees of freedom. 

Since in general the interaction kernel K and the full quark propagators S F are unknown quantities we make two 
(formally covariant) approximations: 

• The propagators are assumed to be of the free form Sf = i{ jl — TOj + ie ) _1 with effective constituent quark 
masses rrij that we use as free parameters in our model. 

• It is assumed that the interaction kernel only depends on the components of p and p' perpendicular to P, i.e. 
K(P,p,p') — V(p_i_p,p±p) with p±p := p — (p ■ P/P )P (instantaneous approximation). 

Integrating in the bound state rest frame over the time component p° and introducing the equal-time (or Salpeter) 
amplitude 



P=(M,0) 



(3) 



we end up with the Salpeter equation |l3( ], which constitutes the basic equation of our model: 

dy v(f,pWY 



$(p) =+Ar(p)7° 

-A+(p) 7 ° 



(2?r) 3 M + CJl + L0 2 

d 3 P ' v&p'W) 

(2?r) 3 M - u)i - uj 2 



7°A 2 -(-p) 



(4) 



Here Af(p) = (wj ± 7°(7P + m i ))/2u>i are projectors on positive and negative energy solutions of the Dirac equation 
and w, = sjp 1 + rnf denotes the kinetic energy of the quarks. 

The simultaneous calculation of the meson masses M and the Salpeter amplitudes $ results by solving the corre- 
sponding eigenvalue problem of eq. (0) with an adequate potential ansatz (see [pi for details). 



III. MODEL INTERACTIONS AND MESON MASS SPECTRA 



The global structure of the experimental mass spectrum reflects a linearly rising confinement force which not only 
produces the Regge trajectories AT 2 oc J but also the energy mass gaps between the radial excitation states. Further- 
more, the confinement force should be flavour symmetric because one finds for every isovector state an energetically 
degenerate isoscalar partner in the experimental mass spectrum. The best known example is of course /o(770) and 
w(782), but also /ii(1170) and &i(1235) and many other pairs, up to ae(2450) and /e(2510). The pseudoscalar mesons 
7T, r], rj exhibit a mass splitting which is of course not compatible with this rule; therefore one has to introduce 
a flavour dependent residual interaction in order to get an appropriate description for these mesons. We will use 
the instanton induced effective quark interaction discovered by 't Hooft |]l6| which in the present formulation acts 
for mesons with vanishing total angular momentum (J = 0) only. In part A of this section, we will focus on the 
confinement force alone and compare the resulting bound state masses of model A and model B to the corresponding 
experimental resonances with J > 0. Presenting the nonrelativistic reduction of both models in part B and discussing 
spin-orbit effects in part C, the residual 't Hooft interaction will be added in part D in order to cover the pseudoscalar 
and scalar mesons. 



A. Confinement Potential and Mesons with J > 



In this subsection, we will discuss two different versions for the confinement force. These two models differ by their 
Dirac structures f®I whereas the coordinate space dependence is chosen to be linear in both parameterizations: 
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'' P ':Vc(p,p')$(p')= I ^V c ((p-p') 2 )mp')r . (5) 



(2tt) 3 ' * ' J (2tt) 3 

Here Vc((p — p') 2 ) is the Fourier transform of the linearly rising potential V(\x q — x q \) = a c + b c ■ \x q — x q \; the 
confinement offset a c and its slope b c are free parameters of our model. So, for given Dirac structure and constituent 
quark masses in a physically reasonable range (m„ w 300 — 400 MeV, m s w 500 — 600 MeV) the calculation of the 
complete J > mass spectrum only depends on these two parameters. In order to fix them, we perform a fit to the 
experimental Regge trajectories and, after some fine tuning of the quark masses, we end up with the parameters of 
model A and model B shown in table ffl.Fl 



1. Dirac Structure F ® F = |(1 ® 1 - 7 <g> 7 ) (Model A) 

This combination of a scalar and a timelike vector Dirac structure has already been used in earlier works ||,[l0 19 
where a very good description of the Regge trajectories M 2 oc J was achieved. The complete (J > 0) spectrum up 
to J = 6 with all its radial excitations up to 2.5 GeV is shown in figs. (JxJ) — (j^) where none of the new experimental 
resonances mentioned in the introduction were used in the fit but only the masses listed by the Particle Data 
Group [^tJ. Consequently a quantitative statement concerning the quality of the higher excitation calculations is 
not possible on the basis of these data alone. Especially, the higher radial excitations of the isovector states 61, a\ or 
of their isoscalar partners hi, fi are not contained in [E7J. Fortunately, the situation has changed illustrated by the 
data shown in table [n] and table [II . Here we have listed several new resonances which all have been found during 



the last few years by various groups and collaborations |l|-0 . 

A comparison of the newly observed resonances with our calculations shows that the Dirac structure 1 ® 1 — 7 <E> 7 
(model A) produces masses that are roughly 150 - 350 MeV too high compared to the experimental values. The 
relativized quark model of Godfrey and Isgur jl5| provides a similar tendency although their deviations are smaller. 
However, this does not exclude a quarkonium interpretation of these states around 1700 MeV and 2100 MeV. Starting 
with the ai(1640) or ai(1700) Q as the well established 2P nh state, a qq classification is much more natural 
than other conceivable interpretations of these resonances: Then the mass centroid of the 2P multiplcts is around 
1700 MeV and the multiplet partners of the ai are expected nearby in mass provided that splittings due to spin-orbit 
and tensor forces are small even for orbital angular momentum L = 1 (P-wave) and L — 3 (F-wave). Indeed, the 
Crystal Barrel collaboration observed a new state, the a2(1660), in the reaction pp — > Tpr\r\ at 1.94 GeV/c f| . A 
similar resonance has been reported in j^] where the authors performed a combined K-matrix analysis of the GAMS, 
Crystal Barrel and BNL data. In addition there is also experimental evidence for the next higher multiplet: The 
authors in 0] studied the process pp — > /2(1270)7r in the mass range 1960-2410 MeV and they observed a 1 ++ and 
a 2 ++ resonance around 2100 MeV which can be identified with the ai(2100) and the 02 (2100), respectively. They 
also found evidence for a 1++ resonance at 2340 MeV, (ai(2340)), a 4++ resonance at 2260 MeV, (04 (2260)), and 
two 3++ resonances at 2070 MeV and 2310 MeV, (<23(20 70) and 03(2310)). However, as already mentioned, neither 
model A nor the model of Godfrey and Isgur produces masses around 1700 MeV and/or 2100 MeV which fit to the 
above quantum numbers. In more detail, model A predicts the first radial excitation of the (pure P-wave) ai at 1876 
MeV and the second radial excitation at 2374 MeV. For the first radial excitation of the 02, the model provides the 
dominantly F-wave state at 1879 MeV below the dominantly P-wave state at 1931 MeV. 



In the isoscalar sector (see table III), the classification of calculated and experimental resonances is more difficult due 
to the doubling of states by the additional ss pair. However, the approximate flavour symmetry of the inter-quark 
forces suggests to expect the first nonstrange radial excitation masses of the /1 and the /2 also around 1700 MeV 
and the corresponding second radial excitations around 2100 MeV. The third nonstrange radial excitation of the fi 
should appear in a mass range of about 2300-2400 MeV and the first nonstrange radial excitation of the /4 around 
2300 MeV. Moreover the isoscalar ground state partner of the 03 should appear around 2000 MeV and its first radial 
excitation partner around 2300 MeV. In fact, the authors in jij found evidence for resonances with these quantum 
numbers. They studied the process pp — > n n°r] for beam momenta of 600 - 1940 MeV/c, corresponding to center 
of mass energies 1960-2410 MeV, and found a 1 ++ resonance at 2340 MeV, the /i(2340) and a 4 ++ resonance at 
2320 MeV, the / 4 (2320). The masses of the isoscalar partners of the a 3 were stated at 2000 MeV, (/ 3 (2000)), and at 
2280 MeV, (/a(2280)). As in the isovector sector, we fix the mass centroid of the 2P nn multiplet around 1700 MeV 



1 Thc parameters of the 't Hooft interaction are fixed in the pseudoscalar sector and will be discussed in part D of this section. 
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assuming that the /i(1700) is the isoscalar partner of the ai(1700). Although the /i(1700) was below the accessible 
mass range in 0, the assumption above is supported by the mass degeneracy of the observed /i(2340) with its 
isovector partner ai(2340). Now, due to the flavour independence of the confinement interaction used in our model, 
the calculated isoscalars are ideally mixed and the nonstrange part of the spectrum is degenerate in mass with the 
corresponding isovector spectrum. Therefore, also the isoscalar masses calculated in model A are too high compared 
to the experimental values. 

In order to fix the centroid of the 2P ss multiplct, we consider the following well established splittings of hn and ss 
isoscalar ground states listed by the Particle Data Group @: [a; (782), $(1020)], [^(1170), fri(1380)], [772 (1645), 
?72(1870)], [a;3(1670), $3(1850)]. These mass splittings of approximately 200-300 MeV leads to an expected centroid 
of the 2P ss multiplet at roughly 1950-2050 MeV provided that the 2P nn multiplet appears around 1700 MeV as 
discussed above. Furthermore the Particle Data Group [[i7| states two isoscalar 2 ++ resonances, the /2(1950) and 
/2(2010), which fit to this energy region such that an arrangement of one of these states into the 2P ss multiplet seems 
to be natural. However, the mass calculations in model A for the strange part of the isoscalar 2 ++ sector provide a 
dominantly F-wave fi at 2148 MeV and a dominantly P-wave ji at 2165 MeV as can be seen in fig. ([j]). Moreover, 
the model produces the first strange radial excitation state of the /1 at 2129 MeV. So, the strange part of the isoscalar 
excitation spectrum comes out too high in model A just like the nonstrange excitation spectrum discussed above. 

In view of these facts, a quarkonium interpretation for the higher mass resonances seems still plausible but obviously 
can only be achieved with a new confinement force. 



2. Dirac Structure r ® T = \{1 ® 1 - 7 5 <g> 7 5 - 7^ ® 7m ) (Model B) 



From first principles QCD hardly gives any clue on the Dirac structure of confinement. Consequently a phenomeno- 
logical study of different Dirac structures is needed. The force with I ® 1 - 7 ® 7 has been found by purely 
phcnomenological arguments. Here the aim was to reproduce the global structure of the meson mass spectrum, 
especially the well established Regge trajectories. Now we focus on some symmetry properties of the strong interac- 
tion. For example the spontaneous breakdown of the approximate chiral symmetry of light flavour QCD leads to an 
effective theory, the well known chiral perturbation theory, which interprets the eight lightest pseudoscalar mesons 
as Goldstone bosons and fixes their interactions by symmetry requirements (apart from phenomenological coupling 
constants). Furthermore it is believed that, due to the absence of a ninth light pseudoscalar meson, the axial U(l)- 
symmetry is explicitly broken by instantons. As we implement this by using the instanton induced 't Hooft interaction 
(see part D of this section) it is interesting to assume a confinement force with axial U(\)-invariance in the present 
framework. We choose the structure 

r $ r = 1(1® I -7 s ® 7 5 -7^07^) . (6) 

Due to the L0i(l)-invariance, this structure obviously leads to parity doublets in the calculated meson spectrum if 
one neglects the quark mass terms in the Salpeter equation (see fig. (Q)). 

The above Fierz invariant combination of a scalar, pseudoscalar and vector part was also investigated by Bohm et 
al. (see j55|) as well as by Gross and Milana (see j56|) but, unlike the investigations presented here, without any 
connection to the complete experimental meson spectrum. 

As figs. (0)-(||) show, this Dirac structure yields ground state Regge trajectories M 2 oc J that are as good as the 
corresponding results achieved with the scalar timelike vector force discussed above. However, both models differ 
essentially in their results for the higher mass resonances: In contrast to model A, the calculated higher excitation 
masses in model B agree remarkably well with the newly observed resonances (see tables || and [II). Here, deviations 



are in general less than 100 MeV. For instance, the first radial excitation of the a\ appears at 1718 MeV (see table [n|), 
which means a lowering of about 160 MeV compared to the corresponding mass calculated in model A (1876 MeV). 
The second radial ai-excitation appears at 2099 MeV, so it is even 275 MeV lighter than the corresponding mass in 
model A (2374 MeV), thus showing excellent agreement with the observed resonance ai(2100) 0. For the third radial 
ai-excitation model B provides 2412 MeV compared to 2791 MeV in model A and 2340 ± 40 MeV the corresponding 
value from experiment Q]. As can be seen in table |l| not only the ai-excitation masses are lowered in model B but 
also all other higher isovector resonances listed there. Due to the flavour independence of the confinement force, the 



same effect occurs in the complete nonstrange and strange isoscalar spectrum (see tabic [II). 

In connection with this general lowering of excitation masses in model B, we can confirm a recent observation found 
by A. V. Anisovich, V. V. Anisovich and A. V. Sarantsev M concerning the systematics of qq states with respect 
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to their radial excitation number. The authors found that for given quantum numbers many mesons fit to linear 
trajectories as 

M 2 = Af 2 + (n - l)^ 2 (7) 

with Mo the mass of the ground state meson [n = 1), n = 2,3,... numerating the radial excitations and u 2 the 
trajectory slope parameter whose value is suggested to be in the region [P 1 = 1.25 ±0.15 GeV 2 . Figs. (||) and (0) show 
several resonance positions taken from || and our corresponding calculations of model A and model B. As expected 
from the discussion of table || and table III, the trajectory slope \j? as given in Q is much too flat to parameterize 



the calculations of model A whereas the model B calculations fit for many mesons to the formula in eq. (|7j) with slope 
fi 2 = 1.25 ± 0.15 GeV 2 remarkably well.0 

In summary, the calculated excitation masses of model B are in general roughly 150-350 MeV lighter than the corre- 
sponding masses of model A, hence in a much better agreement with the newly discovered experimental resonances. 
As we will show below, this significant mass lowering is related to the strong coupling between positive and negative 
energy components of the Salpeter amplitude in model B. 



B. Nonrelativistic Reduction of the Salpeter Equation 



In this subsection, we will show that the negative energy Salpeter components play an essential role for the description 
of higher excitation states in the meson mass spectrum. It is known from earlier works (see e.g. P,|lO|]) that these 
components are very important in calculations of electroweak observables, especially for deeply bound states. There 
it has been shown that a fully relativistic treatment in the Salpeter framework improves the description of many 
observables (for example: pseudoscalar decay constants, two photon decay widths) by orders of magnitude. In order 
to define the negative energy components of the Salpeter amplitude <&, we look at the nonrelativistic reduction of the 
Salpeter equation. This reduction of eq. (^) is reached by neglecting all momentum dependent terms in the energy 
projectors (formally by performing here the limit m t oo) and by expanding the kinetic energies u>i with respect to 
\p i \/m l < 1: 



hm Af(p) = hm = -(1±7 ) =: V 

m— *oo TYh%— *oo ZuJi 2 

2 rrii rrii 



(8) 
(9) 



Inserting these approximations, the Salpeter equation (4) becomes: 

d 3 P > v(p,p')Hp') 



Hp) 



(2tt) 3 M + oj 1 + u> 2 
d 3 p' V(p,p')<f>(p') 



(2tt) 3 M - Qi - Co 2 
Now, if we write the 4 x 4-matrices $ and V§ in block matrix form as 

[P) - { $21 (p) <S>22(P) 



7°P- 



(10) 



2 However, a more detailed study of our calculations with respect to eq. (fit) would prefer especially for the 7r(10 + )- and 
a (10 ++ )-trajectories (fig. (§)) a somewhat larger slope, namely /i 2 (tt(10" + )) = 1.59 GeV 2 and ^ 2 (a (10 ++ )) = 1.56 GeV 2 . 
Then the experimental resonances n (1300 ± 100) and tt(1800 ± 40) and the a (2000+i° ) would fit much better to their 
corresponding trajectories. Consequently, due to the higher slopes, the parameterization in eq. (Q) would lead to higher mass 
predictions, namely tt(2189), tt(2526) and a (1588), a (2375) instead of tt(2070), tt(2380) and ao(1520), a (2260) as stated in 
M. The higher masses then would correspond to the model B calculations 7r(2195), 7r(2496) and ao(1665), ao(2395). 
A universal trajectory fit to our mass calculations provides an average slope of fi 2 = 1.42 ± 0.27 GeV 2 which is roughly 14% 
larger than the value fi 2 — 1.25 ± 0.15 GeV 2 suggested in Q. However, this deviation should not be overinterpreted because 
there is no obvious reason to demand an approximately unique slope for all trajectories. 
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and 



v(p p'mv 1 ) - ( <y*)*®?) <v*)mp') \ 

v ( p,p mp ) - y (yfyn&p) (v<s>) 22 (p,p') ) > 



where each component $,y and (V<$>)ij, i,j S {1,2} is a 2x2 matrix and apply the projectors V in Dirac represen- 
tationf] on both sides of eq. (|Io|), we get: 



$12 (p) 

$ 2 i (p) 



r d 3 p' (V$)2i(p,p") 
J (2tt) 3 M-fiSi+iSa 



f rfV (V$)i2(p,p") 

J (27T) 3 M 



For weakly bound states with M « mi + m 2 one has 

1 



« Ar K - (11) 



M + UJ\ + U>2 M — Cb\ — U>2 



so that the component $21 can be dropped. The above matrix equation then decouples with respect to $12 and $21 
due to the smallness of $21- Therefore one can interpret <I> ++ := $12 as the upper (or positive energy) component of 
the Salpeter amplitude <I>. Consequently, $ := $21 can be interpreted as the lower (or negative energy) component 
of $.0 One should emphasize that the above approximate decoupling for weakly bound state solutions is a consequence 
of the Salpeter equation itself and independently from the special form of the potential operator V in Dirac spacef]. 
With usual constituent quark masses of 300-400 MeV, one therefore can expect to find different Dirac structures 
which all produce for example the p-meson mass around 700-800 MeV with the same potential parameters. However, 
for deeply bound states (M <C TOi + m 2 ) on the one hand and higher excitation states (M 3> TOi + m 2 ) on the other 
hand, neglecting the negative energy components $ is no longer justified. So, in that case even in the nonrelativistic 
reduction the Salpeter equation does not a priori decouple with respect to <I> ++ and $ (see fig. (g)). Since distinct 
potential Dirac structures act on these components in different ways they should also differ in the calculation of 
especially higher excitation masses. This effect then can be used in order to find an appropriate Dirac structure for 
the confinement potential. 



1. Nonrelativistic Reduction of Model A 



The special Dirac structure of model A does not mix positive ( < & ++ ) and negative ($ ) energy components of the 
Salpeter amplitude the components $ ±: F even vanish completely as can be seen by explicit calculation: 

r$r = i(i$i- 7 °$ 7 °) = ( $ _ . 

Therefore the only coupling between $++ and $ is produced by the structure of the Salpeter equation itself due to 
the off-diagonal terms 7°7P in the energy projectors A ± . In the nonrelativistic reduction, these terms vanish and the 



3 Dirac representation of the 7 matrices: 



1 0\ ^ / a\ 5 (01 



7 ^ -1 J ' 7 [ -S j ' 7 I I! !) 

4 The special projector structure of the Salpeter equation (4) allows to express the diagonal components <E> +_ := <E>n and 
$ _+ := $22 in terms of $ ++ and $ as 

$ + " = +Cl $ ++ S-C2 s (12) 

$" + = -Cl $""S + C2 s$ ++ , (13) 

with the shorthand notations s = ap, Ci = ^ im ^ UJ2mi ' ) ( see [@ f° r details). 
5 This statement is at least valid for the Dirac structures discussed in the present paper. 
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Salpeter equation with this Dirac structure then decouples with respect to <J>++ and $ irrespective of calculating 
weakly bound states (M w m-i+m^), deeply bound states (M <C mi+m2) or higher excitation states (M ^> mi+n^): 



<S>++(p) 
(p) 





rfV Vc((p-p') 2 ) 

(27r) 3 (-M)-iii-(Zi 2 



rf 3 P ' Vc((P-P') 2 ) <& + + (p') 

(27r) 3 M— £Di— cDa 



or equivalently 



(H$ ++ )(p) = +A/$++(p) 



(14) 
(15) 



where the operator Ti is defined by: 



(16) 



For physical reasons, the potential parameters of V c should guarantee the positive definiteness of the operator Ti, 
i.e the eigenvalues M of Ti are all positive. Therefore the negative energy component $ has to vanish and we 
end up with the equation for the positive energy component $ ++ which is nothing else but the usual nonrelativistic 
Schrodinger equation for two spin * particles moving in a spin independent confinement potential. 



2. Nonrelativistic Reduction of Model B 



Whereas the structure i(l ® 1 — 7 <8> 7 ) decouples with respect to 3>++ and $ , the combination |(1 <£> 1 — 7 s ® 
7 5 — 7^ ® 7^) leads to relative strong couplings: 

t^t, i/ir^^r 5^5 \ / $++\ 1 \ 1 / -CT^ + CT ct\ 

r$r = -(i$i-7 5 <i>7 5 -7^7,)=( v(i) - J - 2 ( $++ J + 2 U* ++ - -^ + -- J • 

Here we have combined the timelike vector part of 7 M (8 7 M with the scalar part such that the first term on the right 
hand side is nothing else but the result from the Dirac structure of model A. The second term, which arises from the 
pseudoscalar part, swaps the positions of $++ and <f> and of <I>^ and <f> h , respectively. The same exchange of 
components is generated by the vector part but with an additional left-right multiplication by the Pauli-matrices a. 

In the nonrelativistic reduction the components <f>^ and $ h vanish and for each spin (S = 0, 1) the Salpeter 
equation can be written as a system of two coupled 2x2 matrix equations (see Appendix): 

[W(* ±:t + $ T=F )] (p) = ±M$ ±± {p) for 5 = (17) 
[W($ ±=b - <S> TT )] (p) = ±M$ ±± (p) for 5=1 (18) 

where the operator 7i is defined as in eq. (|i"6|).p| 



6 Eq. (|17|) and eq. (y_8[) exhibit an interesting symmetry concerning spin singlet and spin triplet solutions: For given mass ±M 
on the right hand side, a spin singlet solution ( < 1 >++ , "J" -- ) can be transformed into a spin triplet solution with the same mass 
by the transformation 

($++,$—) — ► (±$ ++ ,=f$— ). (19) 

So, in the nonrelativistic reduction the very special combination of coefficients in the Dirac structure 1 — 7 5 ® 7 5 — 7 M ® 7m) 

leads to a mass degeneracy of singlet and triplet spin states although there is an explicit spin dependent term (7 ® 7) in the 
potential. As an example, fig. (JsJ) shows the (spin triplet) p-mass and the (spin singlet) 7r-mass as a function of the parameter 
rh := tuning the nonrelativistic reduction of the Salpeter equation by increasing the quark mass in the energy projectors 
A* indicated by m^. 
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The matrix equations eq. (|17]) and eq. (|1§|) show that even in the nonrelativistic reduction the Dirac structure 
produces nonvanishing negative energy components $ . This can be also seen in fig. ( |l0|) where the positive and 
negative energy components of the ground state a\ radial amplitudes are shown. The same calculation was done with 
the structure i(l eg) 1 — 7 (8)7°) (fig- @) and as expected here the negative energy component completely vanishes 
in the nonrelativistic reduction. The fully relativistic calculation of course produces nonvanishing negative energy 
components also in model A. Whereas the magnitudes of the ground state a\ amplitudes in model A and model 
B roughly coincide, the corresponding amplitudes of the first radial ai-excitation differ substanti ally in both models 
(see fig. (|IT|), fig. ([lj)). Especially the negative energy component produced in model B (fig. ([14)) is much larger 
than its corresponding counterpart in model A (fig. (|ll])). This strong coupling between negative and positive energy 
components in model B has of course an effect on the calculated bound state masses. An instructive illustration of 
this effect is shown in fig. (p^): Starting with the kinetic: energy we have separated the total mass of the first and 
second radial ai-excitation into the potential expectation values with respect to the positive energy components only 
and with respect to all components. As expected from the magnitudes of the corresponding radial amplitudes the 
mass lowering caused by the negative energy components is much larger in model B than in model A. In case of the 
second radial a\ -excitation this effect results in a mass difference of about 300 MeV. 



C. Spin-Orbit Effects 



Finally, we want to point out the most outstanding feature of the Dirac structure |(1 <£> 1 — 7 5 <£> 7 5 — 7 M eg) 7 M ), 
i.e. the generation of additional spin-orbit mass splittings due to the nonvanishing components $ ±= F of the Salpeter 
amplitude $. As we will see in the next subsection, these additional spin-orbit effects in model B provide a completely 
different interpretation of the scalar meson ground state nonet compared to the interpretation proposed by model A. 
In order to study the spin-orbit effects, we decompose the Dirac structure |(1 eg I — 7 s ® 7 s — 7^ eg) 7^) with respect 
to the components <j> ±= F and introduce the spin-orbit parameter a £ [0, i] as follows: 

r$r( a ) = -(i$i-7°$7°)-«(7 5 *7 5 -7$7)= I $ - )-«($++$+-)+«( -<&+-<?)■ 



If a = 0, the components <J> ±=F vanish and the Dirac structure of model A is reproduced, i.e. T $ T(0) = i(lcgl — 7°cg> 
7 ); for a = |, the above expression coincides with the structure of model B, i.e. T $ T(i) = |(Jeg>I— 7 5 ®7 5 — 7 M <8>7 At ). 
In figs. ( |l4| ) and (|l5|), the ground state masses of the ai(l ++ ), ao(0 ++ ) and the A"i(l + ), Aq(0 + ) are shown as a 
function of the spin-orbit parameter a. For vanishing a, no spin-orbit effects are observed whereas for a — i the 
resulting a\-ao and Ai-ATg mass splittingF] add up to 280 MeV and 220 MeV, respectively. This splitting then yields 
an isovector(isoscalar) ground state a (f ) at 944 MeV and an isodoublet ground state K£ at 1100 MeV. So, due 
to additional spin-orbit forces, model B produces the basic scalar ground state mass centroid at roughly 1 GeV in 
contrast to roughly 1.3 GeV in model A. In what follows we will see that this lowering of about 300 MeV combined 
with the instanton induced 't Hooft interaction leads to a completely different nonet classification of the scalar mesons 
in both models. 



D. 't Hooft's Interaction and Mesons with J — 



As it stands the confinement force discussed above is responsible for the global structure of the experimental mass 
spectrum. Especially the well established Regge trajectories can be described by this force alone. However, the 
comparably large mass splittings in the pseudoscalar sector require an additional residual interaction which should 
feature an explicit flavour dependence in order to describe the ^-ry'-mixing and the rj-n mass difference of about 400 
MeV. In fact, 't Hooft and others computed a flavour dependent effective quark interaction from instanton effects 
p6| , p6|j27| whose point-like two body part can be written as |l7j 



7 As the calculation was done without the 't Hooft interaction, the a-dependence of the /i(l ++ ) and /o(0 ++ ) completely 
coincides with the behaviour of their isovectorial partners ai(l ++ ) and ao(0 ++ ). 
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A£(2) = — 22 X, 9eff(i)eikie imn [{^k,c k 1 tfn.cj (*i,c, 1 *m,c m ) + (*fc,c* 7 5 *n,c„) (*J,c, 7 5 *m,cj] 



j k,l c k' c l 



(20) 



where i, k, l,m,n € {u, d, s} are flavour and c&, q, c m , c„ £ {r, <?, 6} colour indices. The e-tensors explicitly show that 
this force only acts on antisymmetric flavour states. Furthermore, due to the positive sign in the Dirac structure 
(1 (g) 1 + 7 5 ® 7 5 ), the [/,A(l)-invariance is explicitly broken such that the [/^(lj-problem is in principle solved by this 
interaction. Note that due to the point-like nature and specific Dirac structure, the Instanton Induced Interaction 
(III) in the above formulation acts on states with total angular momentum J = only. The lowest order contribution 
of this interaction to the Bethe-Salpeter kernel can be extracted as p7fl 



(2tt) 3 



V III (p,p')$(j>')=4G(g,g') 



TZ A (p,p') (Itr [*(p ')] + 7 5 tr [Hp'h 5 ] ) 



(21) 



where 1Z\ represents a regularization function and G{g,g') is a flavour matrix while summation over flavour indices 
is understood. The coupling strengths g (acting on a nonstrange qq pair) , g' (acting on a nonstrange-strange qq pair) 
and the finite effective range A = Am are treated as free parameters in our model. We fix them in the pseudoscalar 
sector in order to reproduce the experimental n-K-rj-r)' mass splitting. The resulting values for the 't Hooft parameters 
in model A and model B can be found in table Q and the corresponding pseudoscalar mass spectra up to 2.5 GeV 
are shown in left part of fig. (|l(|). For a more detailed discussion of the pseudoscalar mesons and their electroweak 
properties we refer to pl| ; here we only want to comment on the r;(1295) and 77(1440): both, model A and model B 
produce only one qq state in the mass region 1200-1500 MeV of the isoscalar '"-spectrum as can be seen in fig. (|l6|). 
The calculated masses (1533 MeV in model ^4, 1446 MeV in model B) would prefer an identification of the 77(1440) as 
the first radial excitation of the ground state 77 such that the 77(1295) is out of place in our quark model. Indeed, there 
are strong experimental hints for a non qq interpretation of this resonance [t30|]5l[ p| due to its absence in the reaction 
pp — > 7r + 7r~7r + 7r~77. Furthermore, disregarding the 77(1295) in the discussion of the 77,77' radial excitation spectrum, 
the bulk of observed and predicted isoscalar '"-resonances fits to the trajectory-like behaviour M 2 = M$ + (n — l)/i 2 
in model B, with an average slope /j 2 (00~ + ) = 1.73 ± 0.08 GeV 2 (see the right hand side of fig. (17)). 



Now, for the sake of completeness concerning the nonrelativistic reduction of our model, we also present this limit for 
the 't Hooft interaction: 



$++(t5) 
$—(jp) 



V 



AG j 







d 3 P ' TC A (f7,p")tr[$++(p , )+$ 



A n r d 3 P ' K A (p,p')tr[$ 
4 G J (2^p" _ 



M- 



M+Ui +CJ2 







So, even in the nonrelativistic reduction, the 't Hooft interaction couples positive and negative energy components as 
can be seen from this expression. Furthermore the term proportional to the identity in eq. (|2l]) vanishes such that in 
the nonrelativistic reduction the 't Hooft interaction does not act for scalar but only for pseudoscalar mesons. In fact, 
the same instanton induced interaction has been used before in a nonrelativistic description within the framework 
of the Schrodinger equation^, both for mesons and baryons 14 1; a satisfying description of the low lying hadronic 
mass spectrum, especially the splitting of the pseudoscalar nonet, was obtained. However, the quoted nonrelativistic 
model is not able to provide a proper description of the scalar nonet due to the absence of any singlet-octet mixing 
mechanism in this sector. 



8 The author in [ p0|j5l| ] suggest the 77(1440) to be the first radial excitation of the ground state 77 compatible with our model 
calculations. Furthermore, in a forthcoming paper we will show, that model B predicts partial decay widths of the ?7(1440) 
into ao(980)7r, K* K and r\a in excellent agreement with the experimental values found in [5l| . 

9 In the framework of the Schrodinger equation, the negative energy components $ were completely neglected. 
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1. Scalar Mesons in Model A 



Now, from an earlier work [ |19fl we know that the present relativistic framework exhibits a natural way to cure this 
shortcoming. There it has been shown that the light scalar mesons can be interpreted as conventional qq states, with 
a small SU(3) mixing angle, governed dynamically by 't Hooft's instanton induced interaction. The calculations were 
done in confinement model A such that the mass centroid of the scalar meson nonet was found at roughly 1.3 GeV 
as discussed in the previous subsection. Then the 't Hooft interaction, fixed in the pseudoscalar sector, provided an 
almost SU(3) octet at about 1400 MeV and a low lying ST/ (3) singlet at 1000 MeV. In more detail, fig. |[|) shows 
the generated mass splitting of the scalar flavour nonet due to the 't Hooft interaction in model A. The dominantly 
singlet states are lowered whereas states with dominantly flavour octet structure are pushed to higher masses by this 
interaction. The calculated ground state masses are predicted to be ao(1321), Tfp(1426), /o(984) and /o(1468), so 
they imply a qq interpretation of the following experimental resonances (see |4j): a (1450), TQ(1430), / (980) and 
/o(1500). A detailed discussion of this interpretation is given in p9| . 



2. Scalar Mesons in Model B 



In confinement model B the corresponding states appear at masses which are roughly 200-300 MeV lighter due to the 
additional spin-orbit effects in this model (see fig. (|l9|) ) : 

a (1057), X *(1187), /o(665), / (1262). 

Before we will compare the masses of these states not only with the resonances positions of jl7| but also with the 
K-matrix pole analysis of V. V. Anisovich and others (see for instance [^2[ ), we make some general remarks concerning 
the scalar nonet classification: 



• General remarks: 

— First of all one should emphasize that, apart from the higher radial excitation states discussed in the 
previous subsection, nowhere else in the calculated mass spectrum the differences between model A and 
model B are as large as in the scalar sector discussed here. This is a very interesting observation, because 
the scalar sector is the only area in the experimental spectrum where mass differences of about 300 MeV 
do not necessarily imply a dismission of one of these models. The reason is the abundant number of 
experimental scalar resonances combined with many different interpretations of their physical nature (a 
detailed and instructive overview is given in J3(]|n| and also in |32]). Especially the scalar isoscalar 
sector features a highly complex resonance structure which to explain is still the most challenging task in 
meson spectroscopy. For instance, the Particle Data Group jlTj lists the following four scalar isoscalar 
resonances up to 1500 MeV: / (400 - 1200) (or cr-meson), / (980 ± 10), / (1370±^g), / (1500 ± 10). In 
order to understand the internal structure of these resonances, an attempt to interpret them as members 
of the l 3 Po and 2 3 Po light-quark nonet is of course the first step one has to do. For that, some kind of 
mass centroid has to be given around which the members of these nonets then split up due to different 
quark flavours, that is to say, due to different quark mass contents and, of course in the present model, 
due to different influences of the instanton induced 't Hooft interaction on these states. As the isodoublet 
mesons A"g(1430) and ^(1950) are the least controversial of the light experimental scalar mesons and as 
they contain one strange and one nonstrange quark, it seems to be obvious to use their masses as fixpoints 
for the l 3 Po and 2 3 Po quark-antiquark nonet, respectively, independent from any theoretical description. 
This scale then favours ground state nonet interpretations in the mass region of about 1200-1600 MeV 
which is in roughly agreement with the energy scale in model -4-Fj 



10 For example, Amsler and Close ji2| built a reasonable scalar nonet in the quoted mass region with the ao(1450) and Kq (1430) 
setting the mass centroid, and their widths setting the scale of the nonet widths. 

Assuming the scalar ground state nonet above 1.3 GeV, Lee and Weingarten |53|j5^ ] found in their lattice gauge theory 
calculations the /o(1710) to be composed mainly of the lightest scalar glueball. They claim that ++ -resonances below 1.3 GeV 
are irrelevant to glueball spectroscopy. 
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— However, V. V. Anisovich [p2| claims that the quoted mass region is noticeably higher than the average 
masses of other mesons which are candidates for the corresponding scalar nonet. In more detail, the 
experimentally well established resonances a (980) and /o(980) are often excluded from a discussion of the 
scalar meson ground state nonet due to their small masses compared to the above scale. In addition, their 
decay properties are often entitled as unusual with respect to 'ordinary' qq mesons which then lead to mainly 
non-qq interpretations of these resonances^. In this context the designation 'ordinary' qq mesons means 
pure nonstrange ++ 3 Po qq states which could be a natural first assignment for these resonances. For 
instance, Godfrey and Isgur ||l5|| argue that this pure nonstrange classification has at least two unpleasant 
consequences in their 'relativized' one-gluon-exchange-plus-linear-confinement potential model: Firstly, 
their model does not account for the experimental fact that both states have sizeable couplings to KK 
final states, suggesting a large strange quark component. Secondly, their model provides strong decay 
widths r(/ (980) -> tttt) rs 400 - 1000 MeV and r(a (980) -> t]tt) k 500 MeV which even overestimate 
the experimental total widths |47) r*°* p (/ (980)) = 40 - 100 MeV and r*°' p (a (980)) = 50 - 100 MeV 
by at least one order of magnitude. In addition, their model predicts two photon decay widths (3^] for 
/o(980) — > 77 and ao(980) — > 77 which are 5-8 times larger than the corresponding measured values. As 
the quoted model describes the decay properties of the vector and tensor qq states with satisfying success, 
the /o(980) and ao(980) were suggested as appropriate candidates for non-qq interpretations. 



Comparison with resonance positions of ^ 

— In contrast, the present fully relativistic quark model with instanton induced forces allows a qq interpreta- 
tion of either the /o(980) (in model A: the /o(984); see |l^] for a detailed interpretation of this state) or 
the ao(980) (in model B: the ao(1057)), i.e. we can not account for a qq interpretation of both states in a 
single confinement model. The model calculations of the total strong decay widths (the detailed treatment 
will be presented in a forthcoming paper @)0 provide r tot (/ (980)) = 126 MeV and r* ot (a (980)) = 55 
MeV in good agreement with experiment. Furthermore, the invariant coupling ratio r(a (980)) = g^ K /g^ n 
is predicted to be r = 1.21 which excellently fits to the latest K-matrix pole analysis of Crystal Barrel 
and CERN-Munich data (r = 1.05 — 1.35) done by A. V. Sarantsev ^o|. Last but not least, we find 
r(ao(980) — > 77) = 0.50 keV EJ in model B, which is in reasonable agreement with the experimental 
estimate T exp (ao(980) — > 77) = (0.30 ± 0.10) keV given in ]4^] . In summary, the model B calculations do 
support a qq interpretation of the ao(980); even more, the results in this model may support its identification 
as the isovector member of the basic l 3 Po qq multiplet. 

— In the isoscalar sector, model B yields a very low-lying dominantly singlet state at 665 MeV which may be 



identified with the broad structure Jo (400 — 1200) (or a- meson) 47 and a state with dominantly flavour 
octet structure at 1262 MeV suggesting an identification with the observed /o(13701^yq) [EtA. Furthermore, 



11 For instance, Weinstein and Isgur argue that the ao(980) and /o(980) can be understood as KK molecules j3^,^8[. The mass 
degeneracy and their proximity to the KK threshold seems to require that the nature of both states must be the same. Also T. 
Barnes |3!| sides with these authors; he claims that the KK molecule picture of the ag (980) and /o(980) is supported by their 
small two photon decay widths. F. E. Close, Yu. L. Dokshitzer, V. N. Gribov et al. |46| suggest to interpret the ao(980) and 
the /o(980) as new types of vacuum excitations ('vacuum scalars') which correspond to quark-antiquark pair creations below 
the Fermi surface. The Jiilich group |^l[ shows that scattering and production data on the ao(980) and /o(980) can be fitted 
by a sum of s-channel and i-channel exchanges without the need for genuine resonances at the KK threshold. On the other 
hand, Morgan and Pennington find a /o(980) pole structure characteristic for a genuine resonance of the constituents and not 
of a weakly bound system that is to say, the extremely attractive 1 = KK interaction may not support a loosely 

bound state, whereas the 1=1 KK interaction is weak and may generate a KK molecule, the ao(980). 
However, there are also attempts to interpret both the /o(980) and the ao(980) as members of the qq nonet: Tornqvist 
claims that these resonances have very large virtual components of KK in their wave functions. However, in order to fit the 
available data on the ao(980), /o(980), /o(1300) and Kq(14:30) mesons as a distorted ++ qq nonet, he has to use a lot of 
parameters (5-6) plus an ad hoc introduced form factor which simulates the overlapping wave functions in the corresponding 
hadronic decay processes. Minkowski and Ochs identify the states /o(980) and ao(980) together with the /o(1500) and 
Kq (1430) as the members of the scalar ground state nonet. They claim that this assignment is supported by phase shift 
analyses of elastic and inelastic -kit scattering as well as recent analyses of pp annihilation near threshold. 

12 In |37]], we will present the interference of two different strong decay mechanisms: instanton induced decay contributions 
from six quark interactions which were studied in [^| for the first time, and quark loop contributions which not only occur for 
scalar and pseudoscalar mesons but also for mesons with J > 0. 
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also model B (as does model A) accounts for the /o(1500) |Q as a qq statej^ and not a glueballQ as can 
be seen in the right part of fig. ([l6]) where the complete scalar excitation spectrum up to 2.5 GeV is shown. 
Whereas in model A the calculated mass at 1468 MeV is identified with the dominantly octet state of the 
scalar ground state nonet, model B provides a dominantly singlet state at 1554 MeV which then belongs 
to the first excited scalar nonet. However, in view of the complexity in this sector that arises from strong 
decay channel couplings and possible mixtures with gluonic or other exotic states, one has to regard also 
decay observables for a more realistic interpretation of the scalar mesons. Some results concerning the 77 
decays of these states can be found in [ [i"l| ; a more detailed discussion which also includes numerical results 
on strong decay widths will be given in a forthcoming paper (see p7|). 

— The classification of the calculated scalar isodoublet states with respect to the listed i0 + -resonances of 
the Particle Data Group |47]] faces a problem in model B: compared to the PDG value of the lightest 
scalar kaon, Mk* ~ 1430 MeV (see ^J), the corresponding model ground state appears at 1187 MeV. 
As discussed above, this lowering of about 250 MeV is a direct consequence from the additional spin-orbit 



forces in confinement model B (see fig. (15)) 



Comparison with K-matrix poles from Anisovich et al.: 

— V. V. Anisovich and coworkers |22 24 4^ |45|] suggested to identify the qq states not with mean resonance 
positions but rather with the poles of the K-matrix fits to the relevant data sets, the so-called 'bare 
states'. They emphasize that pure quark model calculations (as for instance the present one) do not take 
into account the resonance decay, that is to say, these calculations neglect any effects of decay-channel 
couplings on the meson masses. Therefore there are good reasons to compare our calculated masses to 
bare states. Unfortunately, up to now such a K-matrix analysis has not been performed for all quantum 
numbers, so in many meson sectors the bare states are unknown such that an overall comparison with 
our masses is not possible. Moreover, from our point of view in many sectors significant differences in 
mass shifts between bare states and real resonances would be unpleasant due to the fact that our model 
calculations agree well with the global structure of the experimental mass spectrum. 

— However, stimulated by the problems mentioned above with respect to the scalar nonet classification just 
in this sector a multitude of experimental data sets has been analysed by K-matrix pole techniques. The 
most interesting result of these analyses is, that the members of the bare scalar ground state nonet do not 
appear around 1.3-1.4 GeV but in a mass region which is roughly 200-300 MeV lighter than this scale, 
namely 



a£ are (960±30), K* bare {l200tT 5Q ), /o are (720 ± 100), / t,are (1260±^°). 

In particular, the lightest scalar bare kaon appears 200 MeV lower than the amplitude pole. And in fact, 
the calculated scalar ground state nonet of model B shows a remarkable agreement with these bare states 
as can be seen in fig. (p^|)f^] 

— In addition to the scalar ground state nonet, the scalar first excited nonet of confinement model B 

a (1665), X *(1788), / (1554), / (1870) 



13 Again we refer to the paper ]37j for a more solid explanation of a qq interpretation of the /o(1500). 

14 The interpretation of the /p (15001 as the ground state glueball mixed with close-by conventional scalar mesons is strongly 
favoured by Amsler and Close |}l2| . 

1 The low absolute value for the scalar mass centroid of about 1.0-1.1 GeV in model B is not really astonishing due to the 
L0i(l)-mvariance of the Dirac structure |(1 ® 1 — 7 s j 5 — 7^ g) 7 M ) : As discussed above, in the limit of vanishing constituent 
quark masses this invariance leads to parity doublets in the meson mass spectrum. So, in this model the mass splitting between 
states which only differ in their parity quantum number is mainly caused by the nonvanishing quark mass terms in the Salpeter 
equation. The ground states of the strange vector mesons K* and K\ then provide an appropriate estimate of this mass 
splitting. Their calculated masses differ by approximately AM m 400 MeV which also fits to the splitting between the ground 
state masses of the nonstrange parity partners p and ai. Finally, neglecting the Ua ( l)-violating 't Hooft interaction the mass 
of the lightest pseudoscalar kaon K(0~) would appear at roughly 700 MeV (see [|ll|) indicating to expect the ground state 
mass of the parity partner K% (0 1 around 700 MeV + 400 MeV = 1100 MeV, and in fact, this is the value produced by the 
model B calculations (see fig. (|19|)). 
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also coincides with the K-matrix results of the authors in p^- |2^ , ^3| -^5t . They found the following bare 
2 3 P qq nonetQ 

af5 are (1640±40), K* bare (1820tio), / o hare (1600 ± 50), /$° re (1810±ioo). 

Summarizing, the present fully relativistic confining quark model with 't Hooft's instanton induced force as residual 
interaction allows to generate scalar qq states, whose masses and flavour mixings essentially depend on the confinement 
Dirac structure: The c74(l)-violating structure ^(1 <£> 1 — 7 <g> 7 ) fixes the basic mass centroid of the ground state 
nonet at roughly 1.3 GeV, which the 't Hooft interaction then splits into an almost pure flavour sing let /o-state at 
roughly 1 GeV and the flavour octet states /o, do, Kq in the mass region around 1.4 GeV (see |L9|] for a detailed 
interpretation of these states). 

In contrast, the t/A(l)-invariant Dirac structure \ (1 ® 1 — 7 s ® 7 5 — 7^ ® 7 M ) produces the basic mass centroid of the 
scalar ground state nonet at about 1.0-1.1 GeV. Here the 't Hooft interaction lowers the dominantly flavour singlet 
/o-state to roughly 700 MeV and pushes the dominantly flavour octet /o-state to approximately 1.3 GeV. The lightest 
ao-state appears at 1057 MeV and its decay properties are compatible with a qq interpretation of the ao(980). The 
lightest scalar kaon appears roughly 200 MeV lower than the corresponding PDG-resonance ^(1430) j|7|. However, 
as the model presented here does not take into account any effects of decay-channel couplings on the meson masses, 
we followed a suggestion of V. V. Anisovich and coworkers (22| [24|]43| fl5|] , not to identify our calculations with the 
real observed resonances but with the K-matrix poles deduced from appropriate data sets. This assignment then not 
only fits to the calculated members of the ground state scalar nonet but also to their first radial excitation states in 
confinement model B. 



IV. SUMMARY AND CONCLUSION 



Within the framework of a relativistic quark model based on the instantaneous Bethe-Salpeter equation, we have 
studied two different Dirac structures for a linearly rising confinement potential with respect to the complete meson 
spectrum. It was shown that an C/4(l)-hivariant structure of the form |(I(8> 1 — 7 5 ®7 5 — 7^ <8>7 M ) (model B) provides 
an excellent description of the experimentally well established Regge trajectories just as the scalar time-like vector 
combination |(1® 1— 7°<g>7 ) (model A) which has been used in earlier works. However, since the corresponding radial 
excitation spectra were found to be very different in both models we extended the comparison with the experimental 
data to a multitude of new resonances in the mass region 1000-2400 MeV observed in the Crystal Barrel pN- 
annihilation data during the last few years. Whereas the model calculations done with the structure 7°<8>7°) 
overestimate these new data by roughly 150-350 MeV, the combination i(l ® 1 — 7 s <g) 7 5 — 7^ ® 7^) produces masses 
in remarkably good agreement with the newly observed resonances. Here the deviation is in general less than 100 
MeV. Furthermore, the squared masses of the excited states in model B show for given quantum numbers, in contrast 
to model A, a linear dependence on their radial excitation number, M 2 cx n, very similar to the behaviour of many 
experimental resonances recently observed by A. V. Anisovich, V. V. Anisovich and A. V. Sarantsev [0. We thus 
can now relate this observation to the quark-antiquark confinement dynamics in the framework of a constituent 
quark model. We found that the lowering of radial excited states in model B is a direct consequence of the strong 
coupling between positive and negative Salpeter energy components due to the 7 s - and 7 M -part in the structure 
i(l (g> 1 — 7 5 (g) 7 5 — 7 M (8 7 M ). This effect was illustrated by carrying out the nonrelativistic reduction of both Dirac 
structures and it turned out that the Salpeter equation in model A reduces to the usual Schrodinger equation whereas 
in model B the negative energy components do not vanish even in this limit. By virtue of these observations one 
can conclude that any quark model should account for the correct treatment of relativistic effects not only in the 
calculations of deeply bound state masses [M <C m\ + m-i) but also in the case of higher radial excitation states 



There exists a second K-matrix solution for the 2 Po qq nonet: the authors in [|22| [24[[43j [45J find a third scalar isoscalar 
bare state between 1200-1600 MeV, the / 6are (1235 ± 50), which, instead of the / tare (1600 ± 50), can also be used to complete 
the 2 3 Po qq nonet. So, one of the scalar isoscalar bare states in the mass region 1200-1600 MeV is superfluous for the qq 
classification and may be connected to the lightest scalar glueball. The K-matrix solutions lead to positions of the amplitude 
poles in the complex mass plane whose real part may be compared to average mass values taken from the Particle Data 
Group @ ([position of the real part in MeV, PDG 00 Q]): [988±6, o (980)], [1415±25, K*, (1430)], [1015±15, /o(980)], [1300± 
20,/o(1370±f 7 °)],[1499±8,/ (1500)],[1530t^ 
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(M 3> mi + m2). Therefore the mass spectra calculations performed in the fully relativistic Salpeter framework as 
presented in this paper, may be more realistic than the numerous so-called 'relativized' quark model calculations that 
do not incorporate the negative energy components of the constituents in an unrestricted manner. 

A further feature of the fully relativistic Salpeter framework affects the scalar sector: Using the l/^(l)-breaking 
instanton induced 't Hooft interaction in order to compute the pseudoscalar mass splittings, the fully relativistic 
framework automatically generates flavour mixings also in the scalar sector caused by the same interaction. Therefore 
the present model allows to compute realistic scalar qq states, in contrast to its nonrelativistic reduction, where the 
't Hooft interaction acts for the pseudoscalar mesons only. It turned out that the masses and flavour mixings of the 
scalar mesons essentially depend on the confinement Dirac structure: In model A the structure |(I ® 1 — 7 ® 7 ) 
fixes the basic mass centroid of the scalar ground state nonet at roughly 1.3 GcV, which the 't Hooft interaction then 
splits into an almost SU(3) singlet at roughly 1 GeV and an almost ST/ (3) octet at about 1.4 GeV. In contrast, the 
c/4(l)-invariant Dirac structure |(1 <S> 1 — 7 5 (g) 7 s — 7 M (g> 7^) in model B produces the basic mass centroid at about 
1.0-1.1 GeV. Here the 't Hooft interaction lowers the dominantly flavour singlet / -statc to roughly 700 MeV and 
pushes the dominantly flavour octet /o-statc to approximately 1.3 GeV. The lightest ao-state appears at roughly 1 
GeV and its decay properties are not in contradiction with a qq interpretation of the ao(980). The lightest scalar 
kaon appears roughly 200 MeV lower than the corresponding experimental resonance in this sector. However, as the 
model presented here does not take into account any effects of decay channel couplings on the meson masses, we 
followed a suggestion of V. V. Anisovich and coworkers, not to identify our calculations with the observed resonance 
positions but with their K-matrix poles ('bare states') deduced from appropriate data sets. This assignment then not 
only fits to the calculated members of the scalar ground state nonet but also to their first radial excitation states in 
confinement model B. However, up to now such a K-matrix analysis has not been performed for all quantum numbers. 
Therefore the 'bare states' are unknown in many meson sectors, such that an overall comparison with our masses is 
not possible so far. Moreover, from our point of view significantly different mass shifts between 'bare states' and real 
resonances would be unpleasant in many sectors due to the fact that our model calculations agree well with the global 
structure of the experimental mass spectrum. A K-matrix analysis in all meson sectors would be very desirable to 
clarify this matter. 
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APPENDIX: NONRELATIVISTIC REDUCTION OF MODEL B 

The action of the Dirac structure ^(1 <g) I — 7 5 (g) 7 5 — 7 M ® 7^) on the Salpeter amplitude 

can be written in the form 

t^t, l / 1 r^ 1 r 5^5 \ ( $++ \ \ 1 / -a$-+a a<$>—a\ 

r$r = -(i$i-7 5 $7 5 -7^7,) = ( $ - o J-2U++ ^-J+n -s**-*) ■ 

For a better understanding of the explicit spin dependence, it is useful to rewrite the left-right multiplication a ® a 
into a spin projector multiplication from the left: 

ff&ff = (3 V s =o - Ps=i) i,j G {+, -} (22) 

where Vs=o an d Vs=\ denote the projectors on spin singlet (S = 0) and spin triplet (S = 1) states with 'Ps=o + '^ , s=i = 
Is where Is is the identity on the spin space. The left-right multiplication 7 <g> 7 then can be written as 



1G 



7*7= (3P S =o- V s =i) ( $++ 



such that we end up with 

1 <= / -2 <J>~+ <I> ++ + <I>~ 

r$r=-(J® J- 7 5 ® 7 5 -7 m ®7m) = 7> s=0 f $ — + $ ++ _ 2 $+- 



^5= 



o $++ - $- 
$—_$++ o 



In the nonrelativistic reduction the components and $ + vanish and for both spins the Salpeter equation can 
be written as a system of two coupled 2x2 matrix equations: 



[n{^ ±±: + $ T=F )] (p) = ±M$ ±± (p) for 5 = 
[W($ ±:t - $ TT )] (p) = ±M$ ±± (p) for 5 = 1 



where the operator H is defined as in eq. (|l6|). 
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FIG. 1. Isoscalar meson spectrum (J > 0). In the middle of each column the experimental resonance positions [47] and 
their errors are indicated by lines and shaded areas; the lines in the left and right part of each column represent the calculated 
masses in model A and in model B, respectively. 
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FIG. 2. Isovector meson spectrum (J > 0). In the middle of each column the experimental resonance positions [47] and their 
errors are indicated by lines and shadowed areas; the lines in the left and right part of each column represent the calculated 
masses in model A and model B, respectively. 
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FIG. 3. Strange meson spectrum (J > 0). In the middle of each column the experimental resonance positions [47] and their 
errors are indicated by lines and shaded areas; the lines in the left and right part of each column represent the calculated 
masses in model A and in model B, respectively. Note that the calculated K\ states are each 2-fold degenerate for spin S = 
and S = 1, indicated by '2'. 
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FIG. 4. The ao-ir ground state mass difference as a function of the constituent quark mass m n calculated in model A (solid 
line) and model B (dashed line). For vanishing constituent quark masses, the ao-n mass degeneracy in model B is a direct 



consequence of the C/A(l)-invariance of the Dirac structure |(1 < 



7" ® 7" 



7 M ® 7m). In model A this degeneracy is not 



observed due to the J7.A(l)-violating scalar part of the structure |(1 ® 1 — 7° <8> 7°). The calculation was done without the 
explicit C/A(l)-breaking 't Hooft interaction. 
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FIG. 5. The p-n ground state mass split ting as a function of the mass ttia := m ■ m n in the energy projectors 
A* = (wa ± 7° (723*+ mA))/2cJA with lua = \J V 1 + m\ calculated in model B. In the nonrelativistic reduction (in — > oo) 
the (spin triplet) p and the (spin singlet) 7r are degenerate in mass due to the special combination of coefficients in the structure 
|(1 ® 1 — 7° ® j° — 7 M ® 7^). The calculation was done without the 't Hooft interaction. 
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FIG. 8. The squared masses of the isovector Regge trajectory members p, p3 and p 5 as a function of the constituent quark 
mass niA in the energy projectors (see eq. (5)) calculated in the nonrelativistic (Schrodinger equation) model [14] (dotted lines), 
in the full Salpeter model with Dirac structure i(!<g>I— 7°<g>7°) (solid lines) and with Dirac structure i(!<g>l — 7 5 (g>7 5 — 7 M ®7 M ) 
(dashed lines). The calculations were done with the single parameter set of the nonrelativistic (Schrodinger equation) model 
[14]. Whereas the nonrelativistic reduction (mA — > oo) of the structure |(1 ® 1 — 7° 7°) coincides with this model (see eq. 
(12)), the negative energy components 3> generated by the structure i(lCg) 1 — 7 5 ®7 5 — 7 M CS>7 M ) do not vanish in this limit 
(see eq. (26) & eq. (27)). They lead to significantly negative mass shifts with respect to the nonrelativistic model calculations. 
The more the meson mass M differs from the weakly bound state condition M w mi + mi (see eq. (8)) the larger are these 
effects. 
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FIG. 11. Positive {TZ + ) and negative (1Z~) energy components of the first excited radial ai-amplitude in the nonrelativistic 
(left) and fully relativistic version (right) of model A. 
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FIG. 12. Positive {TZ + ) and negative (1Z~) energy components of the first excited radial ai-amplitude in the nonrelativistic 
(left) and fully relativistic version (right) of model B. 
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FIG. 13. Expectation values of the first (lower part) and second (upper part) radial ai-excitation in model A (solid lines) 
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a function of the spin-orbit parameter a calculated with 
B = 1300 MeV/fm). For a = 0, the Dirac structure 
^(J ® 1 — 7 ® 7 ) — q(7 5 ® 7 5 — 7 ® 7) reduces to the structure of model A and no spin-orbit splitting is observed; for 
a — 0.5, it coincides with the structure of model B and generates a spin-orbit splitting of about 280 MeV. As the calculation 
was done without the 't Hooft interaction, the isoscalar partners /i(l ++ ) and /o(0 ++ ) show the same behaviour with respect 
to the spin-orbit parameter a. 




FIG. 15. The ground state ifi(l + )- and J^o(0 + ) -mass as a function of the spin-orbit parameter a calculated with the 
parameters of model B (A = -1135 MeV, B = 1300 MeV/fm). For a = 0, the Dirac structure i(l®l-7°(g>7°)-a(7 5 ®7 5 -7(g>7) 
reduces to the structure of model A and no spin-orbit splitting is observed; for a = 0.5, it coincides with the structure of model 
B and generates a spin-orbit splitting of about 220 MeV. The calculation was done without the 't Hooft interaction. 
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FIG. 16. Pseudoscalar and scalar meson spectrum. In the middle of each column the experimental resonance positions [47] 
and their errors are indicated by lines and shaded areas; the lines in the left and right part of each column represent the 
calculated masses in model A and model B, respectively. 
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FIG. 17. The (n, M 2 )-trajectories for the states 77(00"+) and t/(00~ + ). The resonances indicated with the reference number 
[8] are predictions according to the formula M 2 = Mq + (n — l)/i 2 . The resonances without a reference number are listed in 
[47]. On the left hand side the trajectory slopes fi 2 are taken from [8] and the r;( 1295) is involved. However, the situation 
concerning the radial excitations of the 77 and rj is not clear: In our model, we don't have a serious candidate for the r;(1295) 
to be the first radial excitation of the 77. This is compatible with the observations made in [50,51]. Therefore we propose the 
classification of the ij and r\ spectrum shown on the right hand side, where the 77(1295) does not occur. 
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FIG. 18. Schematic splitting ol the scalar flavour nonet with confinement model A (left), with confinement model A and 
instanton induced force (middle) compared to the experimental spectrum interpreted as qq states [47, 28, 29] (right). 
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FIG. 19. Schematic splitting of the scalar flavour nonet with confinement model B (left), with confinement model B and 
instanton induced force (middle) compared to the K-matrix poles stated in [22] (right). 
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TABLES 



TABLE I. The parameters of the confinement force, the 't Hooft interaction and the constituent quark masses in the models 
A and B. 



Parameter Model A Model B 



't Hooft 


9 [GeV 2 ] 


1.73 


1.62 


interaction 


g' [GeV~ 2 ] 


1.54 


1.35 




Am [H 


0.30 


0.42 


Constituent 


m n [MeV] 


306 


380 


quark masses 


m s [MeV] 


503 


550 


Confinement 


a c [MeV] 


-1751 


-1135 


parameters 


b c [MeV/fm] 


2076 


1300 


Spin structure 


r®r 


i(I®I-7°®7°) 


i(I®I-7 S (g>7 B 
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TABLE II. New candidates for isovector resonances and their masses in [MeV]. The first and second columns show the name 
and publication reference of each resonance, respectively; the third column contains the experimental mass and its error; the 
calculated masses in the full Salpeter model (FSM) are shown in the fourth (model B) and in the seventh (model A) column. 
The corresponding results in their nonrelativistic reduction models (NRM) are also shown. The sixth column shows some 
calculations done by Godfrey and Isgur in their relativized quark model [15]. The dominantly P-wave and F-wave states of 
the ei2 are labeled by P and F, respectively. The dominantly F-wave and H-wave states of the cn are labeled by F and H, 
respectively. 



Name 



Ref. 



Exp. Mass 



Model B 
FSM NRM 



Godfrey &: Isgur 



Model A 
FSM NRM 



a (2025) 



[1] 



2025 ± 30 



2071 



2250 



1932 a 



1927 a 



ai(1640) [47] 
ai(1700) [3] 



1640 ± 42 
1700 



1718 



1845 



1820 



1876 



1927 



ai(2100) 
ai(2340) 



[1] 
[1] 



2100 ± 20 
2340 ± 40 



2099 
2412 



2250 
2586 



2374 
2791 



2528 
3073 



a 2 (1660) 

a 2 (1750) 

a 2 (2100) 

a 3 (2070) b 
a 3 (2310) 
a 4 (2260) 



[47] 

[47] 

[4] 

[4] 
[4] 
[4] 



1660 ± 40 



1752 ± 25 



2100] 



Mil 



2070± 20 
2310± 40 
2260 ± 15 



P:1807 
F:1768 



P:2160 
F:2141 

1926 

2247 

F:2341 
H:2315 



P:1845 
F:1986 



P:2250 
F:2351 

1986 

2351 

F:2351 
H:2492 



P:1820 
F:2050 c 



F:2050 c 
2050 



P:1931 
F:1879 



P:2411 
F:2357 

1951 

2401 

F:2451 
H:2402 



P:1927 
F:2002 



P:2528 
F:2584 

2002 

2584 

F:2584 
H:2658 



a As the ao ground state appears at 1321 MeV in model A, we identify this state with the observed ao(1450) [^] and 
not with the low-lying ao(980). Therefore we assign the observed a (2025) |l| to the first radial excitation of the do 
in model A, whereas in model B the stated masses (FSM: 2071 MeV, NRM: 2250 MeV) correspond to the second 
radial excitation. 

b Although the a 3 (2070) seems to be the ground state in the 3 ++ -sector it does not appear in the listings of the 
Particle Data Group. Its mass is much more questionable than the well established ground state masses in all 
other sectors up to J = 6. Therefore this state was excluded from the fit of the model parameters to the experimental 
ground state masses. 

c These masses belong to the same state. Clearly, its identification with the observed 0.2(2100) Q does not account 
for the two underlying Q,2-states, ci2(1660) and 02(1750), listed in the latest PDG-issue p7| . 
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TABLE III. New candidates for isoscalar resonances and their masses in [MeV]. The first and second columns show the name 
and publication reference of each resonance, respectively; the third column contains the experimental mass and its error; the 
calculated masses in the full Salpeter model (FSM) are shown in the fourth (model B) and in the seventh (model A) column. 
The corresponding results in their nonrelativistic reduction model (NRM) are also shown. The sixth column shows some 
calculations done by Godfrey and Isgur in their relativized quark model [15]. The states labeled with nn contain nonstrange 
quarks only; the states labeled with ss contain strange quarks only. The dominantly F-wave and H-wave states of the fi are 
labeled by F and H, respectively. 



Model B Model A 



Name 


Ref. 


Exp. Mass 


FSM 


NRM 


Godfrey & Isgur 


FSM 


NRM 


/i(1700) 


[7] 


1700 


nn:1718 


nn: 1845 


nn:1820 


nn:1876 


nn:1927 








ss:1958 


ss:2021 


ss:2030 


ss:2129 


ss:2012 


/i(2340) 


[7] 


2340 ± 40 


nn:2411 


nn:2586 




nn:2791 


nn:3073 








ss:2681 


ss:2774 


- 


ss:3094 


ss:2954 


772 (2040) 


[7] 


2040 ± 40 


nn: 1997 


nn:2110 




nn:2156 


nn:2266 








ss:2231 


ss:2279 




ss:2424 


ss:2291 


?72(2300) 


[7] 


2300 ± 40 


nn:2318 


nn: 2463 




nn:2593 


ss:2753 








ss:2571 


ss:2641 




ss:2890 


nn:2829 


h (2000) 


[7] 


2000 ± 40 


nn:1926 


nfi:1986 


nn:2050 


nn:1951 


nn:2002 








ss:2128 


ss:2134 


ss:2230 a 


ss:2193 


ss:2074 


/ 3 (2280) 


[7] 


2280 ± 30 


nn: 2247 


nn:2351 




nn:2402 


nn:2584 








ss:2476 


ss:2514 


ss:2230 a 


ss:2684 


ss:2553 


/ 4 (2320) 


[7] 


2320 ± 30 


nn;F:2342 


nn;F:2351 




nn;F:2451 


nn;F:2584 








nn;H:2315 


nn;H:2492 




nn;H:2402 


ss;F:2553 



a These masses belong to the same state. 
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